Excitonic effects in the optical properties of CdSe nanowires 
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Using a first-principle approach beyond density functional theory we calculate the electronic and 
optical properties of small diameter CdSe nanowires. Our results demonstrate how some approxi- 
mations commonly used in bulk systems fail at this nano-scale level and how indispensable it is to 
include crystal local fields and excitonic effects to predict the unique optical properties of nanowires. 
From our results, we then construct a simple model that describes the optical gap as a function of the 
diameter of the wire, that turns out to be in excellent agreement with experiments for intermediate 
and large diameters. 



Nanowires (NW) exhibit a wide range of unique prop- 
erties 111, including tunable band gaps, ballistic trans- 
port U optical anisotropy Q and strong excitonic ef- 
fects [j]. It is therefore not surprising that for the past 
twenty years NW have emerged as one of the most active 
fields of research in material science This growing 
interest is mainly due to their short term promising ap- 
plications 0, 0: HI and is empowered by a strong demand 
from industries for smaller and more effective devices. To 
a large extent, the novel properties emerge from the lat- 
eral confinement of the electrons in the wire, leading to 
the blue shift of the electronic band gap with decreasing 
diameter. 

In this article we focus on the optical properties of 
CdSe NW from a first-principle perspective, based and 
going beyond standard density functional theory (DFT). 
Beside analyzing how optical properties evolve with the 
diameter of the wire, we test some approximations that 
are commonly used for bulk materials and discuss their 
applicability to electronic excitations in ID systems. Sev- 
eral studies of the optical properties of CdSe wires have 
appeared in the past years P-Q. Classical Q and semi- 
empirical [1] methods have been quite successful in de- 
scribing large-diameter wires, but they fail for small and 
medium diameters. Furthermore, none of articles using 
first-principle methods present in the literature d, Q was 
capable of capturing the physics of electronic excitations 
in these confined systems. This was due to the neglect 
of at least one of the following physical effects that are 
fundamental for an accurate description of nano-scale ob- 
jects: 

(a) Crystal local-field effects. When the polarization 
of the light is perpendicular to the long axis there is 
an accumulation of charges at the wire surface, which 
in turn is responsible for the attenuation of the electric 
field inside the NW. This leads to a strong suppression 
of the absorption for light polarized perpendicular to the 
NW axis, increasing dramatically the optical anisotropy 
of the system Q . This huge effect has already been mea- 
sured (3| in polarized photo-luminescence experiments, 
but unfortunately some recent theoretical works 0, [l^ 
still neglect it. 




FIG. 1. (Color online) Imaginary part of the dielectric func- 
tion for the smallest NW {d = 0.73 nm), calculated in the 
RPA with (LF) and without (NLF) crystal local fields, for 
light polarized parallel (-Ey) and perpendicular (-Ex) to the 
wire axis. Note that the absorption is zero for NLF:iJx- 



(b) Electron-electron interaction and excitonic effects. 
These effects are very pronounced in semi-conducting 
NW due to the attenuation of the screening. The strong 
electron-hole binding is responsible for a red shift of the 
fundamental absorption frequency. Due to this effect, 
NW can be seen as cxciton traps, a property that en- 
dows them of great technological interest [J| . 

To study the importance of local-field corrections, we 
calculated the absorption spectra of unpassivated CdSe 
NW of 5 different diameters (0.73, 1.17, 1.59, 1.77, and 
2.01 nm) starting from DFT Kohn-Sham states and ap- 
plying the random-phase approximation (RPA) to ob- 
tain the dielectric tensor (for details on this and all the 
other approximations employed in the following, refer to 
Ref. Illl ). The neglect of local-field effects is equivalent 
to applying Fermi's golden rule, i.e., to treat only inde- 
pendent particle transitions, and completely ignores the 
inhomogeneity in the dielectric response due to the re- 
duced dimensionality of the nano-object. 

All CdSe NW were assumed to be infinitely long with 
periodic length 7.01 A and with their axis parallel to the 
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FIG. 2. (Color online) Dependence of the RPA absorption 
spectrum with the diameter of the NW. For comparison we 
also show the RPA result for the wurtzite CdSe bulk material. 



wurtzite (001) axis. Ground state calculations were per- 
formed using the DFT code ABINIT and the core 
electrons of Cd ([Kr]4dl0) and Se ([Ar]3dl0) were de- 
scribed by Hamann norm conserving pseudo-potentials. 
We chose the Perdew-Burke-Ernzerhof approxima- 
tion to the exchange-correlation functional. Converged 
calculations required a cutoff energy of 20 Ha and a 
1x1x8 Monkhorst-pack sampling of the Brillouin zone. 
Atomic positions were relaxed starting from the bulk 
wurtzite structure. The converged spacing between the 
wires in our supercell approach was at least 7 A. 

The optical spectra of the wire with 0.73 nm, calcu- 
lated at the RPA level including or neglecting local fields 
is shown in Fig. [1] . Without local fields, there is a small 
anisotropy between the absorption perpendicular and 
parallel to the NW axis, comparable to the anisotropy 
in bulk CdSe. However, and as expected, turning on 
local-field effects suppresses completely the low energy 
absorption peaks in the perpendicular direction, render- 
ing the wire almost transparent below 6.5 eV. This optical 
anisotropy, in agreement with experimental results Q, 
clearly decreases with the diameter, but it is known to 
be still relevant for NW with a diameter of 100 nm. 

To observe how the optical properties vary with diam- 
eter, we plot the RPA spectrum (including local fields) 
for five small wires in Fig. [2] . It is expected that the RPA 
gives a poor quantitative result due to the inappropriate 
treatment of electron-electron and electron-hole interac- 
tions, but it is however possible to extract qualitative 
trends thanks to the partial cancellation of these two 
terms (see below). We see that increasing the diameter, 
and thereby decreasing the confinement effect, leads to 
a red-shift of the spectrum, with the absorption thresh- 
old moving towards the RPA bulk value. There is also 
a redistribution of the oscillator strengths, with the first 
peak loosing intensity. 

To go beyond the RPA and to include the relevant 



FIG. 3. (Color online) Imaginary part of the dielectric con- 
stant in the different approximations: RPA, LRC kernel, and 
BSE for the d = 0.73 nm NW. The ALDA results (not shown) 
are almost identical to the RPA. Also shown is the BSE results 
for the d= 1.17nm NW. 



missing contributions, we performed calculations solving 
the Bethc-Salpeter equation (BSE) [l3l and using time- 
dependent DFT [ll| — both in the adiabatic local density 
approximation (ALDA) and with the model long range 
(LRC) kernels ^ derived from the BSE. The BSE ap- 
proach is the state-of-the-art method for calculations of 
optical absorption and gives results in excellent agree- 



ment with experiments [ll|, llj|. For the LRC and BSE 



calculations, one requires as a starting point the quasi- 
particle band structure, usually obtained within the GW 
approximation [lij . Note that the GW approximation 
predicts accurate band gap energies for CdSe, in contrast 
with the systematic underestimation of the gap obtained 
in DFT. Unfortunately, GW calculations arc computa- 
tionally demanding, even for the small wires studied here. 
We therefore used the following approach to obtain the 
quasi-particle corrections: we solved the GW equations 
for our smallest NW {d = 0.73 nm), and then used a sim- 
ple theoretical model [l3| based on crystal field theory to 
interpolate between this result and the GW gap of the 
bulk 



18|. 



The GW calculations were performed with ABINIT, 
applying the standard plasmon-pole approximation. An 
energy cutoff for the dielectric matrix of 9 Ha and the 
technique of Ref. [l^ to reduce the number of unoccupied 
states in the sums over states were employed. Despite 
the use of a cylindrical cutoff for the Coulomb interac- 
tion [i^l, a distance of 10 A between the wires was re- 
quired to converge the GW corrections to the Kohn-Sham 
band structure. These corrections turned out to be quite 
insensitive to the k-point and band index, in agreement 
with other results for Si and Ge NW 2lJ, [221 , and in con- 
trast with findings for graphene nanoribbons (23j . The 
GW gap for the smallest wire was = 4.84 eV, and 

the final form of the model relating the quasi-particle gap 
with the diameter of the wire (d in nm) turned out to be 
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^GWg1p(^) = l-91 + 2-14/d(cV). 

The BSE results (see Fig. obtained using the code 
YAMBO [i^], prove the existence of strong excitonic 
effects, with the excitonic peaks in the visible energy 
range. The excitonic binding energy compensates al- 
most entirely the large blue-shift coming from the quasi- 
particle corrections, and leads to a transfer of the oscilla- 
tor strength from the higher energy absorption peaks to 
the first peak. On the other hand, the LRC kernels de- 
rived from the BSE, that yield results comparable to the 
solution of the BSE for bulk CdSe [l^ , fail dramatically 
for NW — in fact, as shown in Fig. [3] , even the optical 
gap resulting from DFT-(-RPA calculation is in better 
agreement with the BSE result. The attempt to increase 
the excitonic effect by modulating the empirical parame- 
ters of the LRC kernel did not lead to improvements. The 
reason is that the simple LRC approximations are only 
valid for delocalized excitons, but fail to reproduce the 
considerable binding energy of the very localized excitons 
existing in NW. 

The exciton binding energies from BSE calculations, 
for the NW with diameters of 0.73 and 1.17nm are re- 
spectively 1.6 and 1.13 eV. These values are much larger 
than the binding energy in bulk CdSe [1], 5.1 meV, and 
almost twice as large as the ones found in carbon nan- 
otubes [25|. The decrease of the binding energy with 
the NW diameter, also followed by a transfer of the in- 
tensity of the absorption-edge peak to the higher energy 
ones (see Fig. , reflects the weakening of the excitonic 
effects with the increasing of the exciton spatial extent. 
In fact, the exciton radius of bulk CdSe is 5.6 nm A 
nano-object where one of its dimensions is smaller than 
this value will have strongly bound excitons, due to a 
large overlap between electron and hole. Furthermore 
one expects that the larger the confinement the larger is 
the excitonic binding energy. In our BSE results we also 
found spin singlet dark excitons throughout the spectra, 
with the first one appearing at 3.2 eV for the smallest 
wire. This dark exciton may play a role in non-radiative 
decay processes and thus affect the luminescence proper- 
ties. 

Based on our BSE results for the two smallest wires 
and on the experimental optical gap of the bulk [26| . 
we constructed a simple model to interpolate the optical 
gaps of CdSe NW. The form is the same as the previ- 
ous model for the GW gaps, and the formula connecting 
optical gaps and the diameter of the wire turns out to 
be E'^'°if {d in nm) = 1.74 + 1.05/d(eV). We can see, 
from Fig. g] , that this model describes extremely well the 
experimental results [l^ not only for large NW but also 
in the intermediate regime. Finally, from Fig. 0] we con- 
clude that although the RPA results are red-shifted with 
respect to the latter they still conserve approximately the 
correct slope. 

In conclusion, the first-principle calculations reported 
here show how indispensable is the inclusion of local-field 
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FIG. 4. (Color online) Optical gaps calculated with the RPA 
and the BSE compared to experiment [2^. The hnes con- 
necting the RPA and the experimental points are a guide to 
the eye. The green line is the model interpolation between 
the smallest NW and the bulk. 

and excitonic effect to describe quantitatively the opti- 
cal response of CdSe NW. In fact, calculations based on 
the Fermi's golden rule will never be able to catch the 
physics arising from the low dimensionality of the system, 
which is at the heart of the novel properties of NW. We 
have then provided a simple model for the dependency 
of the optical gap with diameter, which is in excellent 
agreement with available experimental results. Further- 
more we have observed the failure of long-range model 
kernels of time-dependent density functional theory for 
these nano-scale systems, in spite of their success in the 
calculation of the optical properties of bulk CdSe. 
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